Unravelling Two-Photon High-Dimensional Entanglement 
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We propose an interferometric method to investigate the non-locality of high-dimensional two- 
photon orbital angular momentum states generated by spontaneous parametric down conversion. We 
incorporate two half-integer spiral phase plates and a variable-reflectivity output beam splitter into 
a Mach-Zehnder interferometer to build an orbital angular momentum analyzer. This setup enables 
testing the non-locality of high-dimensional two-photon states by repeated use of the Clauser-Horne 
inequality. 
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I. INTRODUCTION 

Entangled qubits play a key role in many applications 
of quantum information 1] and quantum cryptography 
2j. An example of a qubit is the polarization state of 
a photon. More generally, a qudit is a quantum sys- 
tem whose state lies in a d-dimensional Hilbert space. 
The higher dimensionality implies a greater potential for 
applications in quantum information processing and this 
explains the continuously growing interest in methods for 
creating entangled qudits. 

Among these methods, spontaneous parametric down- 
conversion (SPDC) appears to be the most reliable one 
for creating entangled photon pairs 0. Recently, sev- 
eral techniques have been used to create entangled qudits 
from down-converted photons. For example, conserva- 
tion of orbital angular momentum (OAM) in SPDC, has 
been used to create entangled states with d = 3 0, Ü| , 
and a timc binning method was employed to realize states 
with d = 11 0- Recently, spatial degrees of freedom in 
SPDC have been exploited to demonstrate entangle- 
ment for the cases d = 4, 8 Q and d = 6 . 

It is well known that useful high-dimensional entan- 
glement can be witnessed by violation of Bell-type in- 
equalities ^(| which also furnish a test of nonlocality for 
a quantum system. However, tests of d-dimensional in- 
equalities for bipartite quantum systems, require the use 
of at least 2d detectors which becomes exceedingly diffi- 
cult (if not impossible) for large d. 

In a previous paper |ll| we proposed an experiment to 
show the entanglement of high-dimensional two-photon 
OAM states, with two detectors only. This scheme in- 
decd allows to verify the existence of high-dimensional 
non-separability, as demonstrated by our subsequent ex- 
perimental results I n -R- e f- E3 we went on to 
use a 2-dimensional Bell inequality to check the non- 
locality of our OAM-entangled photons. In the meantime 
we have realized that this implicitly assumes dichotomic 
variables, a condition that was not fulfilled by the scheme 
proposed in [ïï| . 

In the present paper, we propose an experimen- 
tal scheme to explicitly test the non-locality (namely, 
the useful entanglement) of very-high-dimensional two- 



photon OAM states (d ~ oo), by using just 4 detectors. 
The advantages of our method with respect to those us- 
ing 2d detectors are obvious for d > 2. Additionally, we 
stress that the scheme we propose is designed to realize 
dichotomic observables. The idea is first to project the 
infmite-dimensional two-photon state onto several differ- 
ent four-dimensional subspaces (in order to select dif- 
ferent four-dimensional two-photon states), and then to 
apply the Clauser-Horne (CH) inequality |lj| to each se- 
lected state. It is not obvious a priori whether such 
a scheme will work or not. In fact several legitimate 
qüestions can be raised: (i) Does this dimensional re- 
duction spoil the entanglement of the two-photon state? 
(ii) Do selected four-dimensional states maximally violate 
the CHSH inequality? (iii) Are distinct four-dimensional 
subspaces equivalent? In the rest of this paper we will 
address these qüestions. 



II. THE PROPOSED EXPERIMENT 

Let us describe the scheme of our proposed experiment 
(Fig. 1). A thin nonlinear crystal yields OAM-entangled 
photon pairs, and the two photons (say a and b) are fed 
into two balanced Mach-Zehnder interferometers. Each 
Mach-Zehnder MZ X , (x — a,b) is made of a 50/50 in- 
put beam splitter and a variable-reflectivity output beam 
splitter, indicated in Fig. 1 with BS and VBSa,, respec- 
tively. We denote with t x and r x the transmission and 
reflection coefficients of each VBSa- and assume 



t x = cos6 x , 
r x — i sin 6 X , 



(la) 
(lb) 



where x = a,b and 9 X £ [0, 27r). Such a VBS can be eas- 
ily realized, for example, by exploiting the polarization 
degrees of freedom of the SPDC photons. Type I crystals 
emit photon pairs with a well defined linear polarization. 
Then, the combination of an half-wave plate before the 
Mach-Zehnder and a polarizing beam splitter as output 
BS of the same interferometer, realizes the desired VBS. 
Another possibility is to use use a Fabry-Pérot étalon 
whose mirror separation can be varied, to realize a so- 
called "Lorentzian beam splitter" |Ï4) , which works as a 



2 



VBS. 

In each channel i, (i = 1, 2) of the interferometer MZ a 
(MZt), there is a spiral phase plate SPP (complemen- 
tary spiral phase plate: CSPP), oriented at ai (Pi). In 
the following we shall restrict our attention to the case 
X2 = Xi + 7T) (x = The output channel "1" of the 

interferometer MZ X is coupled to a single-mode fiber F x i 
which sustains the Laguerre-Gaussian mode LGq with 
waist wq. The output ports of the two fibers F a i and F&2 
are coupled with two detectors D al and D;,i , respectively, 
which measure the twin-photon coincidence rate. The 
experimental scheme also comprises two pairs of imag- 
ing systems (not shown in Fig. 1), which image the twin 
photons from the the crystal to the the SPPs, and from 
the SPPs to the input port of the fibers. 




FIG. 1: Schematic of the proposed experimental setup. The 
boxes MZ a and MZ& represent the Mach-Zehnder interferom- 
eters in the path of the photon a and b, respectively. The 
thick grey lines F X i(x — a,b; i = 1,2), represent the single- 
mode optical fibers. Each of them is coupled with a detector 
Dxi- Other details are given in the text. 




vacuum 



FIG. 2: Detailed scheme of the OAM analyzer in the path 
of the photon a. BS denotes a 50/50 beam splitter, VBS a 
a variable-reflectivity beam splitter, ot\ and a.2 represent the 
two SPPs, and Mi,M2,M3 represent three ordinary mirrors. 
The role of M3 is to maintain the input wave function spatially 
invariant after an even number of reflections. Two additional 
azimuthal-independent phases ip a \ and ip a 2 are accounted for. 



A. The spiral phase plates 

A spiral phase plate, shown in Fig. 3, is a transparent 
diclcctric plate with an edge dislocation that can be freely 
rotated around the plate axis [ï^ . Let z be the axis of the 
plate and x the rotation angle. When a light beam with 
transverse profile crosses such a SPP it acquires an 
azimuthal-dependent phase exp(i/(x, </>)) 

V>(x) -» e i/(x ' w V(x), (2) 

where 

gi/(x,« = e iC(4-x) [e i27r£ 9(x - <P) + 6(0 - x)] ■ (3) 

Here x is the two-dimensional position vector x = (x,y) 
in the transverse plane z = const., <p — </>(x) is the az- 
imuthal angle, and C £ M is the phase shift per unit angle. 
In addition, with 6(A") we denoted the Heaviside func- 
tion which is equal to 1 for X > and to zero otherwise. 
Let S(x) be the quantum mechanical operator represent- 
ing the action of a SPP on the arbitrary single-photon 
state 1-0), and let |x) denotes the position-representation 
state of a quasi-monochromatic photon with a given po- 
larization (position state, for short) 

|X> = hj dV-^íqJlO), (4) 

where q = (q x ,q y ) is the transverse photon momentum 
and [à(q), à^(q')] = <5( 2 )(q — q'). It is easy to see that the 
position states {|x)} are orthogonal and form a complete 
basis in the single-photon Hilbert space 

(x|x'> = 5^ (x-x'), (5a) 

Í = j d 2 x|x)(x|. (5b) 

The quantum operator S(x) can be determined in anal- 
ogy with the classical case, by imposing 

(x\§(x)\1>)=eV<x·*H(x), (6) 



n + 




FIG 3: Schematic drawing of a spiral phase plate (SPP) with 
a step index (= phase shift per unit angle) C = h s (n — no)/A, 
where h 3 is the step height, n and no are the refractive indices 
of the SPP and the surrounding medium, respectively, and 
A is the wavelength of the incident light. In this letter we 
assume C = i + 1/2, {£ — 1, 2, . . .}. The orientation angle a 
is indicated. 
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and assuming ^(x) = (x|-0). Then Eq. © can be rewrit- 
ten as 



which implies, together with the arbitrariness of 



(x|5( X ) 



(7) 



(8) 



This equation shows that the SPP operator S(x) is diag- 
onal in the coordinate basis kets and it is unitary since its 
eigenvalues exp(i/(x, <j>)) have modulus 1. If we multiply 
both sides of Eq. (JSJ) by |x'), we obtain 



(x|S( X )|x') 



(x|x') 



e i/(x,0')( x | x '' 



(9) 



where the second line of Eq. ijüj immediately follows from 
the orthogonality of the position states. From Eq. © 
we easily obtain 



S(x)|x> 



e i/(x»| x \ 



(10) 



which shows, together with Eq. ©, that the SPP oper- 
ator S(x) is symmetric. From Eq. I|1UI) it is straightfor- 
ward to write the corresponding transformation law for 
the spatial creation operators ò^(x): 



S(x)ó t (x)5 t ( X ) = ^'^(x). 



(11) 



To conclude this paragraph, we define the complemen- 
tary spiral phase plate, as a SPP that produces a nega- 
tive azimuthal-depcndent phase shift exp(— i/(x, <j>)) on 
a crossing beam: 

V-(x) -> e -i/ ^W. (12) 

Then from the Hermitian-conjugate of Eq. JSJ) it readily 
follows that the CSPP is represented by Sn(x). 



B. The Mach-Zehnder interferometer 

Figure 2) shows a detailed scheme of the Mach-Zehnder 
interferometer MZ a . Photon a enters the Mach-Zehnder 
through the channel "1" and interacts with the first 50/50 
beam splitter. Channel "2" is fed with vacuum. A gen- 
eral two-mode single-photon (TMSP) state at the input 
of MZ a can be therefore written 



d 2 x£(x)àí(x)|0), 



(13) 



where the subscript 1 denotes the channel 1. The beam 
splitter transforms the input TMSP state \ip a ) in a su- 
perposition of TMSP states [ïq. \ip a ) — > \ip a )i'- 



= J d 2 x£(x)[/tàt (x) £| 0)) 



(14) 



where IT/Jui 

Ü^à\(x)Ü = tà\(x,y) ~r*à\(x, -y), 
Ü^à 2 (x)Ü = ra[(x,—y) + t*à\(x,y). 



(15) 



and we assume t = l/\/2, r = i / y/2 for the tr ansmission 
and reflection coefficients of the 50/50 beam splitter, re- 
spectively. We can then write 



d^xiW 
V2 



à\{x,y)\0)+ià\(x,-y)\Q) 



(16) 

The action of the mirror M2 can be described as 
à\(x,—y) — > rà\{x,y), where r — i, which causes \if> a )i ~ * 
\tpa)2, where 



= / d^x 



V2 



à\(x,y)\0)-àUx,y)\0)]. (17) 



After the first beam splitter and the mirror M2, there 
are two SPPs, one per channel, which perform a unitary 
operation on the TMSP state: \ip a )2 — > |V'o)3· Let S(aii) 
the operator representing the SPP (CSPP) in the chan- 
nel i (i — 1,2). Since under S(cti), the operators <z](x) 
transform as 

aj(x) - S(a z )àl( x )SHa t ) = e i /^·«fij(x), (18) 
then the TMSP state after the SPPs can be written as 

£(x) 



\Í>a 



d 2 x. 



V2 



;i/(«i.«at( x )|0) - e ií(Q3 ^)a|(x)|0) 

(19) 

where Eq. has been used. The action of the mirror 
Mi can be described as à\(x,y) — > — r*à\(x, — y), where 
r = i, which causes \ip a )3 ~ > |V , o)4 ! where 



|^a)4 



,2 C(X) 

d x — — 

V2 



ie if(ai '^à\(x,-y)\0) 
-e i/(Q2 ^ } 4(a;,y)|0) 



(20) 



The second beam splitter transforms the creation opera- 
tors according to [Ï3 

à\(x,-y) -> r a à\{x,y) + t* a à\(x,-y), 



à\{x,y) -> t a à\(x,y) -r*à\(x, -y), 



(21) 



where t a = cos9 a , r a = isin6> a . The TMSP state \ip a )s 
at the output of the Mach-Zehnder MZ a can therefore be 
written as 



à\{x,y) 



ï\(x,-y)}\0). 



(22) 
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The role of the last mirror M3 is to maintain the initial 
spatial wave-function £(x) invariant; its action can be 
described as à\(x, —y) — > rà\(x, y), (r = i), which causes 
|^o)5 ~~ > \ipa)6j where, apart an overall phase factor, 



l^«>6 = / d 



"X 



V2 



à\(x,y) 



P a e i /(ai.« +toe i/(a a .*)] àl(x,y)}\0) 



sm t/ a e 



cos y„e 



Í/(«2,0) 



dx 7rl 



( , )S , Vne i/("i,0) _ S i n( 9 ne i/(«2,<A) aj>,tf) 

4(^y)}|o) 

= /" d 2 x£(x) [^(^(x) + A 2 (0)4(x)] |0). 

(23) 

Equation (|23|l shows that for a given x the TMSP state 
\ipa)6 spans, as 9 a varies, a two-dimensional space deter- 
mined by the orthogonal basis {à}(x)|0), a2(x)|0)}. In 
fact, if we define 



A(4>) 



E{<t>) 



1 (é^ 



(24) 



V2 V ei/2 

where fi = f(cíi, 4>), then from Eq. (|23ll it readily follows 



where 



A(cf>) = R(8 a )E(4>), 



R(e a ) = ( cos í a ~ sin fl 6a ), 

v a > \sm6 a cos6»„ I ' 



(25) 



(26) 



is the well known 2x2 rotation matrix. Therefore, as 
9 a varies from to 27r, the state \1pai6 makes a complete 
rotation in the plane {ò{(x)|0), a|(x)|0)}. It is easy to 
see that if we take in account the azimuthal-independent 
phases tp x i, the orthogonal matrix R(9 a ) must be replaced 
by the unitary matrix U(# a ) defined as 



U(0„) 



COS( 



,}<Pal oi. 



(27) 



Now we can repeat for the photon b the very same 
calculation beginning with the state 



d 2 x£(x)6Í(x)|0), 



(28) 



at the input of the Mach-Zehnder MZ& and ending with 
the state \ipb)& at the output of MZ;,: 



|^6>6= / d 2 x^(x) S 1 (^)6Í(x) + J B 2 (^)ò|(x) |0) 



where 

Bi(<f>) = [cos0 a e- i/(,3l ^ ) -sin^ a e- i/(/32 '^]/\/2, 

B 3 (0) = [sin^e-^^ 1 '^ +cos6 a e- if ^ 2 ^)/V2. 

Note that the minus sign in the exponentials is due to 
the fact that CSPPs (instead of SPP) are used in the 
Mach-Zehnder MZ b . 

Since the Ai(<f>) and Bi(<p) does not depend on the ra- 
dial coordinate r = |x|, in the following we shall indicate 
them as Ai{<j>) and (</>), respectively. 



C. The twin-photon state 

The OAM-entangled state of a photon pair emitted by 
a crystal pumped by a LGq làser beam, can be written 



*> oc / d 2 xAp(r)ó t (x)6 t (x)|0), 



(31) 



where A P (r) = LGg(r, wp) describes the transverse pro- 
file of the pump beam, r — |x|, and wp is the beam waist. 
The state |\&) is clearly non-normalizable, therefore we 
use the symbol "oc" instead of "=" . As we shall see, this 
fact does not represent a problem since all the measurable 
probabilities will be properly normalized. Passing from 
the single-mode to the two-modc dcscription introduced 
in the previous paragraph, we rewrite the two-mode two- 
photon (TMTP) state |#) as 



|*) oc / d 2 xA P (r)a t 1 (x)S t 1 (x)|0) 



(32) 

to indicate that both photons enter the channel "1" of 
their respective interferometers. When the photon pair 
crosses both the Mach-Zehnder, undergoes the trans- 
formation |\fr) — ► |*'): 

I*') oc Jd 2 xAp(r)[[A í ( c l>)à\(x) + A 2 (^àl(^)\ |0) 
L(0)6Í(x) + fl a (0)55(x)] |0)} 

]T / d 2 xAp(r)* y -(0)|x) 

ai | bj j 

(33) 

where |x) ai |x)f,j = àj(x)òj(x)|0) denotes a position state 
with the photon a in the channel i and the photon b in 
the channel j, and 



1,2 



%(0) = ÇijAi{4>)Bj{4>), (1,3 = 1,2), 



(34) 



(29) and ç ió =(3-i-j) + i(3i + 3j - 2ij - 4). 
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D. The single-mode fibers 



culate is 



Figure 2) shows that the output channel i, (i — 1,2) 
of each Mach-Zehnder MZ X , (x — a,b) is coupled with 
the single-mode fiber F^ which sustains the Laguerre- 
Gaussian mode LGq with waist wo- For a proper quan- 
tum mechanical description of the fiber we need to intro- 
duce the Laguerre-Gaussian singlc-photon states defined 
as 



I(jJL,V : C) = I (L'T 



\l,p) = j d 2 xLG^(x)|x). 



(35) 



From the orthogonality property of the LG functions, it 
readily follows 



(l,p\l',p') = 



LG^(x) LGj,(x) 



(36) 



SwSpp'. 



When a photon in the arbitrary state |£) is coupled to a 
single-mode fiber, the fiber transforms the input state of 
the photon in the Laguerre-Gaussian state \l = 0,p = 
0) = |0,0) with probability |(0,0|£)| 2 . Since in our 
scheme the output port of each single-mode fiber is cou- 
pled to a detector, the probability Pij(9 a ,9b) that the 
detector D ai fires in coincidence with the detector D^- is 
given by 



P lJ (9 a ,9 b ) <x I d 2 xAp(r)^(0)(O,O|x) M (O,O|x) b , 

2 

d 2 xA P (r)A 2 (r)%(0) , 
since (0, 0|x) O j = (0, 0|x)y = (0,0|x), and 

<0,0|x) = LG°(r, Wo ) = J-L·e- 1 - 3 ^' = A (r). (38) 



(37) 



As the i&ij((f))s do not depend on r, we can factorize 
Pij ($a> &b) by passing to polar coordinates (x, y) — > (r, (/)): 



d 2 xAp(r)A 2 (r)%(</>) 



2tt 



(39) 



drrA P (r)A 2 (r) / 
Jo 

and cast the radial part aside in order to get 

2 



Pij (0a, db) OC 



(40) 



Finally, from Eqs. H3I22I3Q|I we note that, in practice, 
the only elementary azimuthal integral one needs to cal- 



iLfOM-W,")] 



(41) 



= e^-·'^feit - [(1 - e i27r£ )6(Aí - v) 

-(l-e-^e^-M)] (jjl-v)), 

which reduces to the simpler form 

I(p, l + 1/2) = 2ire-^- v W +1 M (l - , 

(42) 

for £ = l + 1/2, where l £ {0, 1,2,.. .}. 



III. THE CLAUSER-HORNE INEQUALITY 

In the previous section we calculated directly the co- 
incidence probabilities Pij(0 a ,9b) from the TMTP state 
I*') at the output of both interferometers. However, to 
proceed further and test the non-locality of the state | \f r ') , 
we have to specify our scenario more precisely. Let us for- 
malize our experiment as follows. There are two parties, 
say Alice and Bob, who share the two-photon entangled 
state 1^) given in Eq. (|31|l . Each one of the two entan- 
gled photons belong to an oo-dimensional Hilbert space, 
namely the OAM Hilbert space. Alice and Bob have two 
distinct measuring apparatuses: M a and respectively. 
Each apparatus M^, (x = a, b) consists of a two-channcl 
Mach-Zehnder interferometer MZ X , with a parameter X 
at the experimenter's disposal, followed by two (one per 
channel) single-mode fibers F^, (i — 1,2). The output 
ports i = 1 , 2 of each are monitored by two detectors 
D x i and T) X 2 respectively. We stress that in this scenario 
the SPPs rotation angles a and (3 are not experimental 
"knobs" that can be changed during an experiment. Dif- 
ferent pairs {a, /3} define different experiments which use 
the same initial two-photon entangled state |\&). In anal- 
ogy with the polarization case, Alice can choose between 
2 different measurements, say A and A' , corresponding 
to two different choices for the varying-beam-splitter "an- 
gles" 9 a and 9' a , respectively. Similarly, Bob can choose 
between B and B' , corresponding to 9b and 9' b , respec- 
tively. Each time Alice and Bob perform a measurement, 
M x (x — a,b) gives the string {x\,X2}, where Xi = 1 
when the detector fires and ce, = when it does 
not. At this point, our scenario is completely defined: 
We have 2 parties (Alice and Bob), 2 measurements (0 X 
and 9' x ) per party, and 2 possible outeomes ({1,0} and 
{0, 1}) per measurement per party. 

Now we can calculate the quantum mechanical pre- 
dictions for the experimental outeomes. These calcula- 
tions were already donc in the previous paragraph, but 
here we want to repeat them in a slightly different way 
in order to display the dichotomic nature of the prob- 
lem. To begin with, we fix for the rest of this paper, 



6 



ai = a, «2 = ol + n and j3i = j3, [3% — (3 + tt. Moreover, 
we fix £ = l + 1/2, where l E {0, 1,2,.. .}. The TMTP 
state |\&') describes the photon pair at the output 
the two interferometers, just before the fibers. As shown 
previously, each fiber projects any input single-photon 
state in the Laguerre-Gaussian state |0, 0). Therefore, 
from Eq. li vil) it readily follows that the two-photon state 
l'S") after the fibers can be written as 

1,2 

|*") cx Ç|0,0) ai |0,0) 6j 

x í d 2 xA P (r)%(^){0,0|x) oi (0,0|x) 6j ·. 

(43) 

As was shown in Eq. I|39fl . h is possible to write 
d x A.p(r)"$>i } (0)(O, 0| X/ ai 

(0, 0|x) bi 



drrA P {r)A 2 {r) x / dç!> (<£) 
o Jo 



(44) 



— Rx Cij(9 a , 0b), 

where the radial integral R does not depend nor on a 
and (3, nor on # a and We can write thcn 



1.2 



(45) 



where iï has been absorbed into the proportionality 
factor and \i,j) is a shorthand for |0, 0) a i|0, Q)bj- At 
this point, it is straightforward to write the normalized 
TMTP state |* o) after the fibers as 



I* 



01], 



1,2 

i, 3 



(46) 



were we have defined the two-photon amplitudes 

r 27T 
<; I 

Ay (9 a , 6b) 



(47) 



c2tx 



The state |\l/oo) is clearly entangled since the coeffi- 
cients Xij are, in general, not factorable. Moreover, 
it belongs to a 4-dimensional Hilbert space, as a two- 
photon polarization-entangled state, since the continu- 
ous variables (r, (f) have been integrated out. In fact, 
all our operations can be summarized in this way: We 



J 



began with an OAM-entangled two-photon state be- 
longing to an infinite-dimensional Hilbert space H^ xo °. 
Then we performed on this state some unitary opera- 
tions which permitted us to span a certain sub-space 
Finally, we projected the transformed state 



0Í H ab 



onto a 4-dimcnsional Hilbert space Tí ab , the two di- 
mensions (per photon) being provided by the two spa- 
tial modes ("arms") of the Mach-Zehnder interferome- 
ter. In this way the entanglement-preserving mapping 
^ooxoo _^ n 2x2 wag accomp ii s hed. We stress that the 
azimuthal integration in Eq. (|47|) clearly shows that the 
final state l^oo) is entangled because the initial state 
jí') from the crystal was entangled, and not because the 
beam splitters in the MZs created the entanglement [2Ï| • 
Now that we have reduced our problem to a 4- 
dimensional one, there are several inequalities at our dis- 
posal to check the non-locality of the state l^oo)- The 
best known are the Bell inequality [2Ï]], the Clauser- 
Horne-Shimony-Holt (CHSH) inequality [22^. and the 
Clauser-Horne (CH) inequality 0. Since we are propos- 
ing an experiment, we choose here to check the CH in- 
equality which, differently from the CHSH inequality, 
does not require the fair sampling hypothesis |2^| to allow 
the use of unnormalized experimental data. In practice, 
an experimenter choose a measurement, say (A, B), and 
repeats it N times (N realizations) obtaining two strings 
{xiki %2k}, (x = a,b; Xik = 0,1) for each realization 
k, (k — 1,...,N). Then, for N 3> 1, the coincidence 
probabilities Pij(9 a , 9b) are well approximated by the co- 
incidence freqüències Fij(9 a ,6b) 



1 N 

Fij(0 a , 0b) = j^Yl Q (^kb jk - 1/2), 



(48) 



k=l 



where is the Heaviside step function. These freqüèn- 
cies are clearly not "absolute" since in a real experiment 
there are always missing outeomes due, for example, to 
detector incfficiencies and to losses. In other words, we 
can say that the experimenter has not access to the nor- 
malized state l^oo); but only to the unnormalized one 
l'S") (|45|l . Therefore, in order write the CH inequality 
in a useful form for an experimenter, we calculate the 
following unnormalized coincidences probabilities 



Pab(0a, 0b) 
Pab {8a, OO) 

P ab (oo,6 b ) 
P ab (oo,oo) 



Pll, 

Pll +P12, 
Pll +P21, 
Pll +P12 +P21 



P22, 



(49a) 
(49b) 
(49c) 
(49d) 



where pij — \Ci 3 (9 a , 9b)\ , and define the Bell-Clauser- 
Hornc parameter S as 



S 



ï) - Pab(0q, O'b) + Pab(0' a , gt) + PgbiKA) ~ PgbK, gg) ~ Pgbi™, 9 b ) 

P ab (<x,oo) 



(50) 
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Then, the CH inequality requires 

S < 0, 



(51) 



for any objective local theory. 

From Eq. (|49al49d(l it is simple to calculate the four 
coincidence probabilities: They are explicitly given in 
Appendix A. They seems complicated but after a careful 
inspection it is easy to see that if we choose a commou 
orientation a — j3 for the SPPs and the CSPPs for the 
two photons, they reduce to the simpler form 



Pn, 



Pab(0O, OO 
Pab(0a,C° 



Pab(0O, OO 

P ab (oo 7 9 b 



Pab(00, OO 



1 2 

- cos^ 
2 

1 

2' 
1 

2' 



(52a) 
(52b) 
(52c) 



With the particular choice of varying-beam-splitter an- 
gles 

0a = 0, 
0'a = T/4, 

b = tt/8, 
01 = 37T/8, 



we achieve the maximum violation S — (\/2 — l)/2 of 
the CH inequality. This result is vàlid for all pairs of 
"external" parameters (a, fi — a). 

This is the main result of this paper. Differently from 
the polarization case, here we have the additional pa- 
rameter a which can be varied from to 27r in order 
to span part of the infinite-dimensional OAM-entangled 
two-photon Hilbert space. Different vàlues of a define 
different experiments and all these experiments give the 
maximum violation of the CH inequality. 

We stress that the condition a = j3 is sumcicnt but 
not necessary to obtain high violation of CH inequality 
in our scheme. In fact, by numerical search, we found 
many pairs a ^ f3 which produces violations bigger than, 
e.g., 0.204. 

IV. DISCUSSION AND CONCLUSIONS 

What is the meaning of the SPP orientation angles 
pair (a, 0)1 In order to answer this question, let us 
summarize our previous results as follows. Consider a 
detection event represented by the string Dij(9 a ,0b) = 
{{ai, 02}, {bi, 62}}- It is not difficult to show that the 
probability Pjj (9 a , 9b) of such an event can be written as 

Pij (0 a , 0b) = | (0, 0| (0, 0\Üi{a, 9 a ) <g> Ú) [fi, 9 b W n ) \\ 

(53) 

where 



ï_1 ' 2 ' (54) 
x = a, 0, 



[Xj = X + (i — l) 71 ": (x — a i 0)} is the operator repre- 
senting the propagation of a photon through the channcl 
"i" of MZ X , and S(xj) i s the quantum-mechanical oper- 
ator representing a SPP oriented at angle Xj H3- From 
Eqs. 143I44|I it follows that when the event D i] (9 a ,9b) 
occurs, the input state \^ m ) is projected onto the state 
\ui(a,0 a ))\üj(j3,6 b )), where 



\ui(a,6 a )) = Uj(a,6 a )\0,0) 
2 

= Y í R lJ (9 a )\S(a J ))/V2, 
3=1 



(55) 



and \S(ctj)) = S^fa^lO, 0). In a similar manner we define 
{üjfàh)) = ^(/3,0 b )|O,O) and |S(/3 2 )> = S[ 
From (S(ai)\S( aj )) = 6 tJ = (S{Pi)\S(fis)) Ü 
lows that {\S(a)), \S(a + vr))} and {\S(f3)), \S((3 + tt))} 
form an orthogonal two-dimensional basis for the pho- 
tons a and b, respectively. Therefore, Eq. (|45(l telis 
us that the state \ui(a,9 a ))\üj(fi,9b)) onto which Al- 
ice and Bob project their state lí' 111 ), is confmed to 
the four-dimensional two-photon subspace spanned by 



the basis {\S(on)) ® \S(Pj)}}, 



1,2). Moreover, 



we can see that, e.g., the basis {\S(a)),\S(a + tt))} 
defines a dichotomic subspace as the basis {|iï),|V)} 
does in the polarization space. It is clear then that 
when we choose a pair (a, (3) of SPPs orientations, we 
uniquely fix a four-dimensional two-photon subspace. 
Since the CH inequalities are applicable to the counts 
at any pair of detectors which measure dichotomic vari- 
ables (irrespective of the dimensionality of the bipar- 
tite quantum state \W n ) under examination [2H 127^). 
we can choose other pairs (a',/3') of SPPs orientations 
(which define other four-dimensional two-photon sub- 
spaces), repeat the measurements, and find again the 
maximum violation of CH inequalities. Now, providing 
that the state vectors {\S( X )), \S{ X + tt)), \S( X ')), \S( X ' + 
tt)) , \S( X ")), | S( X " +7r)>, . . . } ( X = a, (3) are chosen to be 
linearly independent, we can extend the CH test to the N 
pairs {(a,í3), {a',f3'), (a",(3"), («W,^')} defining 
N pairs of two-dimensional subspaces whose union de- 
fine a 2N x 2N two-photon subspace. In this way we can 
demonstrate the non-local nature of the high-dimensional 
two-photon OAM-entangled states. 

In summary, in this paper we proposed a novel ex- 
perimental setup to investigate the non-locality of high- 
dimensional two-photon OAM-entangled states gener- 
ated by SPDC. We use a pair of modified Mach-Zchnder 
interferometers (one per photon), as OAM analyzers. 
Inside each MZ there are two SPPs (one per arm) 
which can rotate around their axes and permit us to ex- 
plore the infinite-dimensional two-photon Hilbert space. 
The output port of each MZ is made of a reflectivity- 
varying beam splitter which acts as a polarizer in the 
two-dimensional space defined by the two spatial modes 
(the two arms) of each MZ. When the output ports 
of these OAM analyzers are fed into single-mode opti- 
cal fibers, the effective dimensionality of the two-photon 
Hilbert space reduces from 00 to 4. Because of this 
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entanglement-preserving dimensional reduction, our ex- 
perimental scheme permits us to check the non-locality 
of the two-photon OAM-entangled state, by using a 
d x N a x Nb = 2x2x2 inequality |2S|. In this way 
we found the maximum violation of the CH inequality 
for any four-dimensional two-photon subspace selected 
by the SPPs orientations. Moreover, because of the 
strict analogy between ours four-dimensional two-photon 
subspaces and four-dimensional two-photon polarization 
space, other interesting experiments (e.g. teleportation 
of spatial degrees of freedom) can be implemented by 
using our scheme. 



P ab (6 a , oc) = 37r 2 + 2,5 2 -^(V + 5| + 2|<S| + |7r-ó|) 
+7T Í2S - |tt + í| + \tt - S\ \ cos(20 Q ), 

(A2) 
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APPENDIX A 



P ab (oo,9 b ) = 37r 2 + 2<5 2 -7r(V + <5| + 2|<5|-r>-<5|) 
+7T 26 - \tt + 6\ + \tt - 6\ cos(2(9 fc ), 

(A3) 



For completeness, we give here explicit expressions for 
the unnormalized probabilities displayed in Eq. I|49al 
I43ÏÏ1 . Here<5 = a-/3. 



sm 2 9b 



2ir\6\cos 2 (0a-e b ) 
2vr 2 + S 2 



cos 



Pab(0a,0b) = S 2 COS 2 (0 a -0b 

+ sin 2 9 a ^n 2 si: 

-2tt (5 + ]tt - 5\) | + cos 2 8 a {-!T 2 cos 2 8 b 
+ sin 2 9 b 2tt 2 + 5 2 - 2tt (-6 + |tt + S\) } 
+- sin(20 a ) sin(20 6 ) [tt (|tt + 5\ + \n- S\) 
-\n + S\\tt - S\ , 

(Al) 



P ab (oo,oo) = 6 7 r 2 +4(5 2 -27r(|7r+ ( 5|+2|(5| + |7r-(5|). (A4) 
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